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We consider double parton distributions in the general case in which the virtualities 
of the interacting partons are different. We elaborate the corresponding evolution 
equations and their extension to next-to-leading logarithmic accuracy. 
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In hadron-hadron collisions it is often assumed that final states containing high-mass systems 
or high transverse momentum jets are generated by a single hard scattering which involves 
one parton from each colliding hadron. The possibility, however, of multiple hard scatterings 

> , 

should be considered as well. One might therefore consider the case in which two hard 

00 

interactions occur within the same hadron-hadron collision as an approximation to the full 

his 

a. 



I. INTRODUCTION 



multiple parton interactions contributions. Several experimental results indeed support t 
possibility and are based on the analysis of the four-jet jl-3| and 7 + 3 jets channels 0, 



Multiple parton interactions has been first modelled and included in modern Monte Carlo 
8|. Very recently detailed phenomenological investigations on double 



^ ! event generators 

parton scattering have appeared in the literature. The y fo cus on the four-jet l9|, [10|, double- 



12 



final states. 



inclusive-forward-pion production [13J], same-sign W and Z plus jets 
The efforts to identify processes which could be maximally sensitive to the contributions of 
double-parton scattering (DPS) is driven by two main interests. On the one hand a carefull 
assessment of phase-space region where DPS events might impact searches for new physics 
is needed. On the other hand a genuine understanding of hadron structure in high energy 
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collisions in terms of multi-partons distributions would emerge from these studies. Most of 
the predictions reported in the phenomenological analysis are based on the simplified model 
in which double parton distributions (DPD) are supposed to be the product of single-parton 
distributions. This assumption is indeed reasonable given the regime of low parton fractional 
momenta presently accessible at hadron colliders. Such an assumption simply disregards any 
longitudinal-momentum and flavour correlation between the two interacting partons from 
each hadron, so that each one evolves according to standard DGLAP equations [ijjj. The 
main virtue of such an approach is that it is technically appealing since numerous single 
parton distributions sets are available. The scale dependence of double-parton distributions 
has been worked out in Ref. [l^. With respect to standard single-parton distributions 
evolution equations (DGLAP), they do contain an additional term which is responsable for 
dynamical correlation between the interacting partons. Quite recently a new set of double 
parton distributions has been obtained by means of numerical integrations of the DPD 
evolutions equations. The initial conditions are such that DPD preserve under evolution a 
number of momentum and flavour sum rules The evolution equations elaborated in 



Ref. [15| however assume that both the interacting parton have the same virtualities. 

Numerical studies Q, [ijl and the arguments given in Ref. [l^ indeed indicate that 
the characterizing scale for double parton scattering is the transverse momentum of the 
final state products. One may therefore consider the production of a gauge boson of mass 
M 2 = Q\ in the first hard scattering associated with jets produced in the second hard 
scattering and characterized by the jet transverse momentum P 2 = Q\. We indicate with 
Qi and Q\ the factorization scales for the two hard processes. The low P 2 regime, with 
P t 2 <C M 2 , for which we expect significant contributions from DPS events, is not covered by 

purpose of this paper is to obtain DPD 
evolutions equations for different virtualities of the interacting partons. Then we consider 
the extension of the formalism beyond the leading logarithmic approximation. By using jet 
calculus rules we work out the inhomogeneous term at next-to-leading order accuracy and 
connect the real two-loops splitting functions arising in DPD evolution equations to the one 
appearing in fracture functions evolution equations at the same level of accuracy. Our main 
results are all framed within the Jet Calculus formalism since it proves to be an efficient 
tool for calculating multi-parton distributions properties and only an ab initio calculation 
could bring these findings on a firmer ground. 
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This paper is organized as follows. In Sec. [TT] we review the basics of Jet Calculus for- 
malism and recover known results on DPD. In Sec. II III we work out the DPD evolution 
equations at different virtualities. In Sec. IIVI we guess the evolution equations for DPD at 
next-to-leading order accuracy. Finally we summarise our results in Sec. |V] 



II. PRELIMINARIES 

The double-parton distributions (DPD) D 3 ^ ,32 (xi,Q 2 ,x 2 ,Q 2 ) are interpreted as the two- 
particle inclusive probability of finding in a target hadron a couple of partons of flavour j\ 
and j2, fractional momenta X\ and x 2 and virtualities up to Q\ and Q 2 , respectively. The 



special case in which Q\ = Q\ = Q 2 has been considered in detail in Ref. 15] . According to 
Jet Calculus [18j . the distributions at the final scales, Q\ and Q 2 ,, are constructed through the 
parton-to-parton functions, E, which themselves obey DGLAP-type evolution equations: 

where Pl(u) are the Altarelli-Parisi splitting functions. Inserting the initial condition 
E 3 (x,Ql,Q 2 ) = S 3 5(l — x) eq. ([T]) can iteratively be solved to give 
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Ei(x, Q 2 , Q 2 ) = 5{6(1 -x) + ^(x) In ^ + 0{a 2 s ) . (2) 
Therefore the functions E provide the resummation of collinear logarithms up to the accu racy 



with which the Pj,{u) are specified. We may therefore express, by Jet Calculus rules 18], 
the double-parton distributions Dj^' 32 (xi,Ql,X2,Q 2 ) as 

Di^( Xl ,Ql,x 2 ,Q 2 2 )= (3) 
r X2 — —Dt J HzuQl,z 2 ,Ql)Ep(^,Ql,Ql)E 3 f(^,Ql,Qi) + 

[ MW{Ql,Ql) d £ f~ X2 d Jl f' Zl ^D^^^^E^^Q^E^^Q 2 ) . 

The first term on r.h.s., usually addressed as the homogeneous term, takes into account the 
uncorrelated evolution of the active partons found at a scale Ql in D 3 ^' 32 up to Ql and Q 2 ,, 
respectively. The second term, the inhomogeneous one, takes into account the probability 
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to find the active partons at Q\ and Q\ as a result of a splitting at a scale /i^, integrated 



over all the intermediate scale at which such splitting may occur. The distribution D°^^ rr is 



a) r>3\,32 



zi + z 2 



r>Ji 



Z\ + z 2 



(4) 



The distributions in eq. (jlj) are the single parton distributions and the ?j' j2 are the real 



Altarelli-Parisi splitting functions 



18]. Both terms in eq. (EJ) are shown in Fig. ([I]). 



Ji 



J 2 

r 2 





FIG. 1: Pictorial representation of both terms on right hand side of eq. (|3|). Black dots symbolize 
the parton-to-parton evolution function, E . 

Due to strong ordering in parton virtualities, the maximum scale in the nl integral is set to 
Min{Q\, Qz)- The scale Ql is in general the (low) scale at which DPD are usually modelled, 
in complete analogy with the single-parton distributions case. In the present context it 
also acts as the factorization scale for the correlated term, since all unresolved splittings for 
which nl < Ql are effectively taken into account in the definition of J2 (zi, Ql, z 2 , Qq). 
The limits on convolutions integrals in eq. (jlj) are fixed by momentum conservation, 



Z\ > X lf Z 2 > X 2 , Zi + z 2 < 1 



(5) 



where Z\ and z 2 are intermediate partons fractional momenta and the last condition guaran- 
tees that their sum never exceeds the incoming hadron fractional momentum. The, lowest- 
order, real Altarelli-Parisi splitting functions P{j 9 (u) and Pjj 9 (u) both contain an infrared 
singularity at the endpoint, u = 1. It is however easy to show that such a singularity is 
always outside the triangle defined by eq. (JSJ) in the [z\,z 2 ] plane, provided that the triv- 
ial condition Xi,x 2 > holds. In the "equal scales" case, Q\ = Q\ = Q 2 , we may take 
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the 
Ref. 

Q 



OKarithmic derivative with respect to Q 2 in eq. and recover the result presented in 

15]: 

a s (Q 2 ) f 1 du 



2 dD^( Xl ,x 2 ,Q 2 



a s (Q 2 
2tt 



dQ 2 
1 du 



2tt 



P^(u)D^ k ( Xl / Ul x 2 ,Q 2 )+ 



u 



-Pf(u)Di^( Xl ,x 2 /u,Q 2 ) + ^(Q 2 )^ (xi+x 2 ,g 2 ) af 1J - a /_x j _ 

2_ U 2,71 Xi + X 2 J \Xi + x 2 



• (6) 



The first and second terms on the right-hand side are obtained through the Q 2 dependence 
contained in the E functions, while the last is obtained from the Q 2 dependent limit in 
the fi 2 integration in the correlated term. The evolution equations therefore resum large 
contributions of the type a s \ia(Q 2 / Ql) and a s ln(Q 2 /fi 2 ) appearing in the uncorrelated and 
correlated term of eq. (j3J), respectively. 

III. EVOLUTION EQUATIONS FOR DIFFERENT VIRTU ALITIES 



Let us now consider the general case in which the partons initiating the two separate hard 
scatterings have different virtualities, Q\ and Q 2 , respectively with Q\ <Q\. The evolution 
equations for the higher scale is obtained by taking the logarithmic derivative of eq. (JHJ) with 
respect to Q\ 



Ql 



dD{^(x 1 ,Q 2 ,x 2 ,Q^ 

1-X 2 



dQl 



1— XI 



1-2 



dffi 



<Xs(fJ? a 



dz\ 

Z l J XI 



' 2 dzi 

XI Z l J X2 

1-21 J~ Eij'f- i _ ..2 



1 ^D^(z u Ql z 2 , Q 2 )Ep (^,Q 2 , Ql) + 

Z 2 J i \ Z\ / 



1 dz 2 F 3 h (z l + Z2 1 til pjjj' a f Zl 

j ' \Zi+Z 2 )^^\zi 



Z'2 



Zi + z 2 
2tt 



ElH-^lQl 



1 -n 2 (u)E^^lQ 2 2 ) 



^2 U 
2 



z 2 u 



(7) 



and using eq. ([!}. Reordering the integrals, we get 
^dD^\x^Qlx 2 ,Q 2 ) <x s {Qf '' 1 



dQl 

dz x r 1 -* 1 <i-- 



2tt 



—Pt(u) 
u 



^D^( Zl , Ql z 2 , Q 2 )El (|, Ql Q\)E\ (^-, Ql Q\) + 
nl - Zl dz 2 F[( Zl + z 2 ^ 2 s) p 3 y 2 



11 d 2 oc s (vi) f u dz i 



Ql 



27T^ 



Z'2 



Zl + z 2 



Zl 



Zi + Z 2 J J i V Zi 



J2 Z 2 U 
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It is now easy to recognize, through direct comparison with eq.OH]), that the term is square 
brackets is the double parton distribution Dj^ ,k (xi,Ql,x 2 /u,Q 2 1 ). The desidered evolution 
equations then becomes 

Q2 dDi^\x^,x 2 ,Ql) = a^ r 1 ^ p n iu)D n,k ixuQlx2/ ^ Qlh (g) 

^ J X2 u 

1 — x 1 

We could obtain the same result in a rather different way. We can in fact exploit the following 
property of the E function 

Ej(x,Ql,Ql)= f ^EH-,Ql,Ql)Ei(u,Ql,Ql) . (10) 

J r£ til V III / 

The latter can be checked, for example, by expanding the E functions in power of a s as 
given in eq. (j2J). By using eq. (fTOI) . eq. © can be recast in the much compact form 

Dt»( Xl , Ql x 2 , Ql) = ^Dt k ( Xl , Ql w 2 , Ql)Ei> Ql Ql) . (11) 



By direct substitution it can be checked that eq. (fTTj) is indeed a solution of eq. (J9J). With 
respect to "equal scale" DPD evolution equations we notice the disappereance of the inho- 
mogenous term. This is due to the fact that the correlations up to a scale Q\ given by the 
inhomogeneous term are taken into account by the "equal scales" evolution equations and 
properly built into D^' k (xi,Ql,w 2 ,Ql) . The evolution of the second parton from Ql to 
Ql is uncorrelated due to strong ordering in virtualities assumed in the leading logarithmic 
approximation. From the numerical point of view therefore DPD at different virtualities can 
be obtained evolving Dj^ ,k (xi,Ql,w 2 ,Qf) with the "equal scale" evolution equations up to 
Qi, eq. dHJ), and then using the latter output as initial condition in eq. for Q\ > Q\. 
We have threfore proven the conjecture put forward in Ref. [16] and actually implemented 
numerically [26] . For completeness we have also considered the DPD evolution equations in 
Q\. Provided that Q\ < Q\ and using the same techniques through which we have derived 
eqs. (jHJ) and (JSJ) we get 

+ aJ M £-»■ ^^^^(^^(^^ . (12) 

In this case the evolution equations contain an inhmogeneous term which arises due to the 
explicit Ql dependence on the integral in eq. ([3]). Since the factorization scale are kept 
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different, the latter does contain explicitely the function E(Q\, Q 2 ), which cannot be further 
simplified. To avoid a direct calculations of the E function, the double-parton distributions 
for unequal final scales should be obtained therefore via the two step procedure mentioned 
above. 



IV. EVOLUTION EQUATIONS TO NLLA 



In this section we address the problem of deriving the structure of DPS evolution equa- 
tions at next-to-leading logarithmic accuracy. The aim therefore is to provide some guidance 
for an eventual ab initio calculation. At present, in fact, such an accuracy is not required 
since, given the scarce experimental information available, we do not even have sufficient 
data to test whether the scale dependence predicted by DPD evolution is supported. Jet 



Calculus techniques has been succesfully extended up next-to-leading 



to improve the perturbative description of time- like parton cascades 20]. For space-like par 



ogarithmic accuracy 



ton cascades instead, which is the case we are actually interested in, the formalism has not 
been extended beyond leading-logarithmic accuracy. However a couple of calculations have 
been performed in the context of semi-inclusive Deep Inelastic Scattering. In particular the 
one-particle inclusive cross sections up to order 0(a 2 ) have been calculated in Refs. 2l|, |22|. 



Such calculations carefully consider hadron production collinear to the hadron remnant 



where the introduction of fracture functions 



23] is shown to be necessary to factorize addi- 



tional collinear singularities appearing in the calculations in that phase-space region. The 
fixed order calculations at 0(a 2 ) allows the authors to derive the fracture functions evolution 
equations to next-to-leading logarithmic accuracy, as well as the two-loop, unknown, real 
splitting functions, P^\ Fracture functions evolution equations can be calculated, as DPD, 



within the Jet Calculus formalism 



24 



25| and they do contain an inhomogenous term as 



well. While in the fracture functions case the partons emitted by the active one hadronizes 
through a fragmentation function, in the DPD one, the emitted parton is allowed to further 
evolve and eventually initiate a second hard scattering. 

When evaluting the evolutions equations at next-to-leading logarithmic accuracy the 
evolution equations for the parton-to-parton functions E must be properly modified to 



1 du 

U L 



P 



(0),*, 



U) + 



a s (Q 2 
2tt 



-P 



(l).i, 



U 



E*(x/u,Qt,Q 2 ), (13) 



8 



where p(°\u) and P^'(u) are the one- and two-loops 19[ Altarelli-Parisi splitting functions, 
respectively. This in turn implies that the two homogenous terms in DPD evolution equa- 
tions in eq. §6§ are modified by adding the two-loop splitting functions contributions. On the 
contrary, the derivation of the inhomogenous term to next-to-leading logarithmic accuracy 
is not trivial so, in the following, we will construct it explicitely in the "equal scales" case. 
The correlated term can be written therefore as 



•^W s 2 )) 



a s (id 



2tT(1 



p(o)i^ 



1 dw n±ir- 

X1+X2 J Xl %1 J X2 

2 



dz 2 

Z-2 



dr x dr 2 dui du 2 



a. u 2 ) ■ ^rr/^iu 



■ Ep^fi^Q 2 ) E^(r 2 ,fi 2 ,Q 2 ) 8{x l -r 1 z 1 ) 5(x 2 -r 2 z 2 ) 5(z x - u t w) 5(z 2 -u 2 w). (14) 

In the above equations PjP^^ 2 (ui,u 2 ) gives the probability that a parton f splits to three 
partons, where the first, j[, and a second, j 2 , have respectively a fraction u\ and u 2 of the 
incoming parton momentum f and the third is integrated over. Integrating the 5-functions, 
which implements longitudinal momentum conservation, one gets 



h,corr 



Xi,X2, Q , Q z 



dw 



xi+x 2 



w 



2 F l ( W ,Vs) 



1—X2 



■*«(*)((!.*.*) + 

J \ W / \ W W / 



dz\ 

Z l JX 2 

Zl z 2 



l-zi 



dz 2 

Z2 



. w w 

- t>: Q 2 ) E% 



t*„ Q ■ (15) 



'i \zi rs ' " J ^ \z 2 

As already noted, the inhomogenous term in DPD evolution equations is due to the explicit 
Q 2 dependence in the upper limit of fi 2 integration. In order to obtain it we set fi 2 = Q 2 in 
eq. flTlj]) . multiply by Q 2 , and use intial condition on E, Ej(x,Q 2 ,Q 2 ) = 515(1 — x). Adding 
the homogeneous contributions, the final result reads 



Q 



2 dD{^(x u x 2 ,Q 2 ) 



+ 



a s (Q 2 
2tt 



<9Q 2 

a s (Q 2 ) 



+ 



xi+x 2 



2tt 

dw 
w 2 



u s (Q 2 

2tt 
du 



du 



*i u 



a s (Q 2 
2tt 



*2 U L 



P^(u) + 



a s (Q 2 
2tt 



-P 



(1)J2, 



u 



D£' k (x 1 ,x 2 /u,Q 2 )+ 



wP$ 



f - >o(k; — Zi — z 2 )H ^ — i' 



(i),iij2/^i ^2 



(16) 
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It should be noted however that the kernels ' 11 2 (w,t>) reported in Refs. [21|, |22[ do 
express the probability that a parton f splits into a parton j[ with a momentum fraction u 
of the incoming parton, into a parton j' 2 with a mometum fraction v of j[, the third being 
integrated over. Therefore they are related to the ones appearing in eq. (fT6"j) by the following 
mapping 



Ml 



Ml, 



Mi 



(17) 



The additional integral in the inhomogeneous term does appear since the momentum is not 
anymore constrained in the 1 — > 2 splitting. The DPS evolution equations to next-to-leading 
logarithmic accuracy for different scales can be obtained by the same arguments given in 
Sec. IIHI We just quote the final result which reads 



Q 



dQl 



2tt 



du 



u 



Pi 0) ' j2 (u) + 



a s{Ql) p (i)j a , 

" J 7, 



2tt 



u 



Dt\x 1 ,Qlx 2 /u,Ql), (18) 



provided that Q\ <Q\. 



V. SUMMARY 



We have considered double parton distributions in the general case in which the two factor- 
ization scales are kept different and derivered the corresponding evolution equations. The 
results of the present calculation support the guess put forward in Ref. 16) and recently 



implemented numerically 



26| widening the range of possible phenomenlogical investigations 



on double-parton scatterings. We have also derived the general structure of the DPD evolu- 
tion equations at next-to-leading logarithmic accuracy and indicated how to transform the 
two-loops real splitting functions present in the literature in order to be used in the present 
context. Both results should be confirmed by performing an ab initio calculation. 
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